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Abstract. 

We describe statistical and coherence properties of the radiation from x-ray free electron lasers (XFEL). 
It is shown that the X-ray PEL radiation before saturation is described with gaussian statistics. Particularly 
important is the case of the optimized X-ray PEL, studied in detail. Applying similarity techniques to the 
results of numerical simulations allowed us to find universal scaling relations for the main characteristics 
of an X-ray PEL operating in the saturation regime: efficiency, coherence time and degree of transverse 
coherence. We find that with an appropriate normalization of these quantities, they are functions of only 
the ratio of the geometrical emittance of the electron beam to the radiation wavelength. Statistical and 
coherence properties of the higher harmonics of the radiation are highlighted as well. 
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1. Introduction 

Single pass free electron laser (FEL) amplifiers startingfrom shot noise in the electron beam 
have been intensively developed during the last decadeq§|. An origin for this development was 
an idea bom in the early eighties to extend the operating wavelength range of FELs to the 
vacuum ultraviolet (VUV) and x-ray bands [1-3]. Significant efforts have been invested into 
the development of high brightness injectors, beam formation systems, linear accelerators, 
and undulators. The result was rapid extension of the wavelength range from infrared to hard 
x-rays [4-13]. The first dedicated user facility FLASH at DESY in Hamburg is in operation 
since 2005 and provides wavelength range from 6.5 nm to 50 nm [14]. LCLS at Stanford 
has been recently commissioned and delivers radiation in the 0.15 nm - 1.5 nm wavelength 
range [13]. The two other dedicated facilities that are under construction at the moment, the 
European XFEL, and SCSS at Spring-8 [15, 16]. 

The high gain FEL amplifier starting from the shot noise in the electron beam is a very 
simple device. It is a system consisting of a relativistic electron beam and an undulator. 
The FEL collective instability in the electron beam produces an exponential growth (along 
the undulator) of the modulation of the electron density on the scale of undulator radiation 
wavelength. The initial seed for the amplification process are fluctuations of the electron beam 
current. Since shot noise in the electron beam is a stochastic process, the radiation produced 
by a SASE FEL possesses stochastic features as well. Its properties are naturally described 
in terms of statistical optics using notions of probability density distribution functions of the 
fields and intensities, correlation functions, notions of coherence time, degree of coherence, 
etc. 

Development of the theoretical description of the coherence properties of the radiation 
from SASE FEL has spanned more than twenty years (see [19-36]. This subject is rather 
complicated, and it is worth mentioning that theoretical predictions agree well with recent 
experimental results [7-10, 37-40]. Some averaged output characteristics of SASE FEL 
in framework of the one-dimensional model have been obtained in [19,20]. An approach 
for time-dependent numerical simulations of SASE FEL has been developed in [21,22]. 
Realization of this approach allowed one to obtain some statistical properties of the radiation 
from a SASE FEL operating in linear and nonlinear regime [23,24]. A comprehensive study 
of the statistical properties of the radiation from the SASE FEL in the framework of the same 
model is presented in [25]. It has been shown that a SASE FEL operating in the linear regime 
is a completely chaotic polarized radiation source described with gaussian statistics. Short- 
pulse effects (for pulse durations comparable with the coherence time) have been studied 
in [22, 26, 27]. An important practical result was prediction of the significant suppression 
of the fluctuations of the radiation intensity after a narrow-band monochromator for the case 

§ Following the terminology of quantum lasers (amplified spontaneous emission, ASE), the term "self amplified 
spontaneous emission (SASE)" in connection with an FEL amplifier, starting from shot noise, started to be used 
in [17]. Note that this essentially quantum terminology does not reflect physical properties of the device. In fact, 
free electron laser belongs to a separate class of vacuum tube devices, and its operation is completely described 
in terms of classical physics (see [18] for more detail). 
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of SASE PEL operation in the saturation regime [26]. Statistical description of the chaotic 
evolution of the radiation from SASE EEL has been presented in [28, 29]. 

The first analytical studies of the problem of transverse coherence relate to the late 
eighties [30, 31]. Later on more detailed studies have been performed [32]. The problem 
of start-up from the shot noise has been studied analytically and numerically for the linear 
stage of amplification using an approach developed in [31]. It has been found that the process 
of formation of transverse coherence is more complicated than that given by a naive physical 
picture of transverse mode selection. Namely, in the case of perfect mode selection the degree 
of transverse coherence is defined by the interdependence of the longitudinal and transverse 
coherence. Comprehensive studies of the evolution of transverse coherence in the linear and 
nonlinear regime of SASE EEL operation have been performed in [33-35]. It has been found 
that the coherence time and the degree of transverse coherence reach maximum values in 
the end of the linear regime. Maximum brilliance of the radiation is achieved in the very 
beginning of the nonlinear regime which is also referred as a saturation point [34]. Output 
power of the SASE EEL grows continuously in the nonlinear regime, while the brilliance 
drops down after passing saturation point. 



2. Operation of an FEL amplifier 



A single-pass EEL amplifier starting from the shot noise in the electron beam seems similar 
to the well known undulator insertion device: in both cases radiation is produced during 
single pass of the electron beam through the undulator. To reveal principal differences, 
we first recall the properties of incoherent radiation. Radiation within the cone of half 
angle ^con = Vl + K'^/ilV^w) has relative spectral bandwidth AA/A ~ 1/-/Vw near the 
resonance wavelength A = (Aw/27^)(1 -I- K'^). Here A^ is the undulator period, is 
the number of undulator periods, 7 is the relativistic factor, K — e\^H^/2T:mc^ is the 
undulator parameter, is the rms undulator field, and m and e are the electron mass 
and charge, respectively. Radiation energy emitted by a single electron in the central cone 
is E, ~ 47r2e2/rM2j/[A(l + K^)]. Here Ajj = 1 and Ajj = [Jo(Q) - Ji(Q)] for a 
helical and a planar undulator, respectively, J„(Q) is a Bessel function of nth order, and 
Q = K'^ /2/{l + K"^). Wavepackets emitted by different electrons are not correlated, thus 
radiated power of the electron bunch with current / is just the radiation energy from a single 
electron multiplied by the electron flux 7/e: 



W^incoh ^ 



A 



1 + 



42 



(1) 



In the free electron laser, electron beam density is modulated by the period of resonance 
wavelength A. Let us consider a model case of an electron beam with the gaussian distribution 
of the current density with rms width (7, and an axial modulation 7(2;) = /[l-|-ainCos(27rz/A)]. 
Total power radiated by a modulated electron beam has been derived in [41]: 



W^coh = 



2c 



1 + 



Aj,N^F{N) , 



(2) 
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Figure 1. Plot of the universal function F (N) given by eq. Q. Dashed curve shows an asymptote for a 
wide electron beam F = \/ {2ttN). 



where 



F{N) 



TT 



1 \ / 4A^2 

arctan ( — — + In 



2N I V4A^^ + 1 



(3) 



= ka'^/L^ is the Fresnel number, k = 2n/X, and = N^X„ is the undulator length. 
In Fig. [H we present the plot of the universal function F{N). It exhibits a simple behavior in 
the limits of large and small values of Fresnel number: F(N) ^1/ {2txN) for ^ oo, and 
F{N) ^ 1 for iV ^ 0. 

Analysis of expressions ^ and (O tells us that incoherent radiation power corresponds 
to the radiation power of the modulated electron beam with effective modulation amplitude of 



«in ~ 1/ V N„IX I (ec) . Note that N^IX / (ec) is the number of electrons on the slippage length 
A^„A. Now we have quantitative answer to the question: how much the FEL radiation power 
exceeds the power of incoherent undulator radiation? In the free electron laser an amplitude 
of the electron beam density modulation reaches values of about unity, and the ratio of the 
radiation powers (coherent to incoherent) is a factor of about the number of electrons per 
slippage length. 

Enhancement of the beam modulation in the free electron laser occurs due to the 
radiation-induced collective instability. When an electron beam traverses an undulator, it 
emits radiation at the resonance wavelength. The electromagnetic wave is always faster than 
the electrons, and a resonant condition occurs such that the radiation slips a distance A relative 
to the electrons after one undulator period. The fields produced by the moving charges in one 
part of the electron bunch react to moving charges in another part of the bunch leading to 
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a growing concentration of particles wherever a small perturbation starts to occur. Electron 
bunches with very small transverse emittance (of about radiation wavelength) and high peak 
current (of about a few kA) are needed for the operation of short wavelength FELs. 

The description of high gain PEL amplification refers to the class of self-consistent 
problems for which the field equations and equations of motion must be solved 
simultaneously. Characteristics of the amplification process can be obtained with a 
combination of analytical techniques and simulations with time dependent EEL simulation 
codes (see [18, 36] and references therein). The amplification process in the SASE EEL is 
triggered by the shot noise in the electron beam, then it passes the stage of exponential growth 
(also called the high gain linear regime), and finally enters saturation stage when the beam 
density modulation approaches unity. In the linear high-gain limit the radiation emitted by the 
electron beam in the undulator can be represented as a set of self-reproduced beam radiation 
modes [42]: 

^ = j ciu;exp[iu;(2:/c - t)] x ^Y^Ankii^, z)^nk{r,uj) e:x.^[Knk{'^)z + incj)] (4) 

n,k 

described by the eigenvalue A„fc(c<j) and the field distribution eigenfunction $„fc(r, uj). Here 
uj = 27rc/A is the frequency of the electromagnetic wave. At a sufficient undulator length 
the fundamental mode (having maximum real part of the eigenvalue) begins to be the main 
contribution to the total radiation power. Erom a practical point of view, it is important to find 
an absolute minimum of the gain length Lg = l/Re(Aoo) corresponding to the optimum 
focusing beta function. In the case of negligible space charge and energy spread effects 
(which is true for XEELs) the solution of the eigenvalue equation for the field gain length 
of the fundamental mode Lg and optimum beta function (3opt are well approximated bjljj]: 



1.67 



/?opt - 11.2 (^) (5) 



/ ; A2/3 KAjj 

1/2 3/2,1/2 

I J \KAjj 

Here is normalized emittance, and I a = mc^/e ~ 17 kA is Alfven's current. 
Dimensionless EEL equations are normalized using the gain parameter V and the efficiency 
parameter p [18]: 



Ia AAw7^ 



1/2 



^ = 17- 

Analysis of the dimensionless EEL equations tells us that the physical parameters describing 
operation of the optimized EEL ©, the diffraction parameter B and the parameter of betatron 
oscillations kjs, are only functions of the parameter e = 27re/A [18,34,44]: 

II General fitting expression including energy spread effects can be found in [43]. 
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B = 2Ta^uj/c^l2.5x e^/^ , 

= 1/{(3T) ~ 0.158 X r=^/2 _ (7) 

Note that Eqs. ^ and ([7]) are accurate in the range 1 < e < 5. 

The diffraction parameter B directly relates to diffraction effects and the formation of 
transverse coherence. If diffraction expansion of the radiation on a scale of the field gain 
length is comparable with the transverse size of the electron beam, we can expect a high 
degree of transverse coherence. For this range of parameters the value of the diffraction 
parameter is small. If diffraction expansion of the radiation is small (which happens at large 
values of the diffraction parameter) then we can expect significant degradation in the degree of 
transverse coherence. This effect occurs simply because different parts of the beam produce 
radiation nearly independently. In terms of the radiation expansion in the eigenmodes ^ this 
range of parameters corresponds to the degeneration of modes [18]. Diffraction parameter for 
an optimized XFEL exhibits strong dependence on the parameter e (see eq. (|7])), and we can 
expect that the degree of transverse coherence should drop rapidly with the increase of the 
parameter e. 



3. Definitions of the statistical properties of radiation 



We describe radiation fields generated by a SASE PEL in terms of statistical optics [46]. 
Longitudinal and transverse coherence are described in terms of correlation functions. The 
first order time correlation function, gi{t, t'), is defined as: 



9i{r,t-t') 



{\E{f,t) \'){\mt') 



1/2 



(8) 



For a stationary random process the time correlation functions are dependent on only one 
variable, t = t — t' . The coherence time is defined as [18,47]: 



\gi{r)Vdr. 



(9) 



The transverse coherence properties of the radiation are described in terms of the 
transverse correlation functions. The first-order transverse correlation function is defined as 



7i(rl,f/x,2,t) 



{E{f^,z,t)E*{Tl^,z,t)) 
(|E(r-l,z,t)|2)(|E(r/^,z,t)|2)''^' 



where E is the slowly varying amplitude of the amplified wave, E = E{f±, z, t) Qi'^oiz/c-t) q 
We consider the model of a stationary random process, meaning that 71 does not depend on 
time. Following ref. [34], we define the degree of transverse coherence as: 

/ |7i(^±, rf±)\'^I{r±)I{ff^)df±dff± 



C 



[//(rx)dr-l] 



(10) 
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where / oc is the radiation intensity. 

An important figure of merit of the radiation source is the degeneracy parameter 5, the 
number of photons per mode (coherent state). Note that when 5^1, classical statistics are 
applicable, while a quantum description of the field is necessary as soon as 5 is comparable 
to (or less than) one. Using the definitions of the coherence time Q and of the degree of 
transverse coherence (fTOl) we define the degeneracy parameter as 



5 = NphrcC, (11) 

where Nph is the photon flux. Peak brilliance of the radiation from an undulator is defined as 
a transversely coherent spectral flux: 

dcu (A/2)' A3 

When deriving right-hand term of the equation we used the fact that the spectrum shape of 
SASE PEL radiation in a high-gain linear regime and near saturation is close to Gaussian 
[18]. In this case the rms spectrum bandwidth o^i and coherence time obey the equation 



4. Probability distributions of the radiation fields and intensities 

The amplification process in the SASE EEL starts from the shot noise in the electron beam, 
then it passes the stage of exponential amplification (high gain linear stage), and finally enters 
saturation stage (see Eig. O. The field gain length of the fundamental radiation mode in the 
high gain linear regime is given by eq. ©, and saturation is achieved at the undulator length 
of about Lsat — 10 X Lg for the parameter space of modem X-ray EELs. Eigures [3] and |4] 
show evolution of temporal and spectral structure of the radiation pulse along the undulator: 
at 0.5Lg (beginning of the undulator), 5Lg (high gain linear regime), and lOLg (saturation 
regime). Eigure [5] shows snapshots of the intensity distributions across a slice of the photon 
pulse. We see that many transverse radiation modes are excited when the electron beam enters 
the undulator. Radiation field generated by SASE EEL consists of wavepackets (spikes [22]) 
which originate from fluctuations of the electron beam density. The typical length of a spike is 
about coherence length. Spectrum of the SASE EEL radiation also exhibits a spiky structure. 
The spectrum width is inversely proportional to the coherence time, and a typical width of a 
spike in a spectrum is inversely proportional to the pulse duration. The amplification process 
selects a narrow band of the radiation, the coherence time increases, and spectrum shrinks. 
Transverse coherence is also improved due to the mode selection process dH). 

Eigure [6] shows the probability distributions of the instantaneous power density / oc 
(plots on top), and of the instantaneous radiation power P oc J I{fi_) d^L (plots on bottom). 
We see that transverse and longitudinal distributions of the radiation intensity exhibit rather 
chaotic behavior. On the other hand, probability distributions of the instantaneous power 
density / and of the instantaneous radiation power P look more elegant and seem to be 
described by simple functions. The origin of this fundamental simplicity relates to the 
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Figure 2. Evolution of main characteristics of SASE EEL along the undulator: brilliance (solid hne), 
radiation power (dash-dotted Une), degree of transverse coherence (dashed Une), and coherence time 
(dotted line). BrilUance and radiation power are normalized to saturation values. Coherence time and 
degree of transverse coherence are normalized to the maximum values. Undulator length is normalized to 
saturation length. The plot has been derived from the parameter set corresponding to e = 1. Calculations 
have been performed with the simulation code FAST [48]. 




t [a.u.] t [a.u.] 



Figure 3. Temporal structure of the radiation pulse at different undulator lengths. Indexes 1, 2, and 3 
correspond to the undulator length of 0.5Lg, SLg, and lOLg, respectively. The plot in the right column 
represents the enlarged fraction of the plot in the left column. Calculations have been performed with the 
simulation code FAST [48]. 

properties of the electron beam. The shot noise in the electron beam has a statistical nature that 
significantly influences characteristics of the output radiation from a SASE PEL. Fluctuations 
of the electron beam current density serve as input signals in a SASE PEL. These fluctuations 
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Figure 4. Spectral structure of the radiation pulse at different undulator length. SoUd, dashed, and dotted 
line correspond to the undulator length of 0.5Lg, 5Lg, and lOLg, respectively. The plot on the left-hand 
side shows only the envelope of the radiation spectrum. The plot on the right-hand side represents an 
enlarged fraction of the radiation spectrum. Calculations have been performed with the simulation code 
FAST [48]. 




Figure 5. Snapshots of the power density distribution in a slice at the undulator length of 0.5Lg (left plot), 
5Lg (middle plot), and lOig (right plot). Calculations have been performed with the simulation code 
FAST [48]. 



always exist in the electron beam due to the effect of shot noise. Initially fluctuations are 
not correlated in space and time, but when the electron beam enters the undulator, beam 
modulation at frequencies close to the resonance frequency of the FEL amplifier initiates the 
process of the amplification of coherent radiation. 

Let us consider microscopic picture of the electron beam current at the entrance of the 
undulator. The electron beam current consists of moving electrons randomly arriving at the 
entrance of the undulator: 

N 
k=l 

where 5(. . .) is delta-function, (-e) is the charge of the electron, is the number of electrons 
in a bunch and tk is the random arrival time of the electron to the undulator entrance. 




Figure 6. Probability density distributions of the instantaneous power density / = \E\'^ (top), and 
of the instantaneous radiation power P (bottom) from a SASE PEL at different stages of ampHfication: 
linear regime, saturation regime, and deep nonlinear regime (undulator length of 5Lg, lOLg, and 15Lg, 
respectively). Solid lines on the power density histograms (top) represent negative exponential distribution 
(fT4l i. Solid lines on power histograms (bottom) represent gamma distribution (fTSl) with M = l/o"p. Here 
e = 2. Calculations have been performed with the simulation code FAST [48]. 



The electron beam current I{t) and its Fourier transform /(a;) are connected by Fourier 
transformations : 



y N 

I{uj) = / e^^*/(t)rft= (-e)^e^^*'= . (13) 

— oc 

It follows from eq. (fT3l) that the Fourier transformation of the input current, I (to), is the sum 
of large number of complex phasors with random phases (pk = ^tk- Thus, harmonics of the 
electron beam current are described with gaussian statistics. 

The FEL process is just an amplification of the initial shot noise in the narrow band 
near the resonance wavelength A when both harmonics of the beam current and radiation are 
growing. An FEL amplifier operating in the linear regime is just a linear filter, and the Fourier 
harmonic of the radiation field is simply proportional to the Fourier harmonic of the electron 
beam current, E{uj) = Ha{oj — ujq)I{uj). Thus, the statistics of the radiation are gaussian - 
the same as of the shot noise in the electron beam. This kind of radiation is usually referred to 
as completely chaotic polarized light, a well known object in the field of statistical optics [46]. 
For instance, the higher order correlation functions (time and spectral) are expressed via the 
first order correlation function. The spectral density of the radiation energy and the first-order 
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time correlation function form a Fourier transform pair (Wiener Khintchine theorem). The 
real and imaginary parts of the slowly varying complex amplitudes of the electric field of the 
electromagnetic wave, E , have a Gaussian distribution. The instantaneous power density, 
/ = l-Ep, fluctuates in accordance with the negative exponential distribution (see Fig. (6]): 



Any integral of the power density, like radiation power P, fluctuates in accordance with the 
gamma distribution: 



wherer(M) is the gamma function with argument M = l/(Tp,and(Tp = ((P — (P))^)/(P)^ 
is the relative dispersion of the radiation power. For completely chaotic polarized light 
parameter M has a clear physical interpretation - it is the number of modes [18]. Thus, 
the relative dispersion of the radiation power directly relates to the coherence properties of 
the SASE FEL operating in the linear regime. The degree of transverse coherence in this case 
can be defined as [18]: 



It is shown in ref. [34] that such a definition for the degree of transverse coherence is 
mathematically equivalent to (flOl) . 

When amplification process enters nonlinear stage and reaches saturation, statistics of 
the radiation significantly deviate from gaussian. Particular signature of this change is 
illustrated in Fig. [6l We see that the probability distribution of the radiation intensity is 
not the negative exponential, and the probability distribution of the radiation power visibly 
deviates from gamma distribution. Up to now there is no analytical description of the 
statistics in the saturation regime, and we refer the reader to the analysis of the results of 
numerical simulations [34]. General feature of the saturation regime is that fluctuations of the 
radiation intensity are significantly suppressed. We also find that the definition of the degree 
of transverse coherence ([T6l) has no physical sense near the saturation point. 

When we trace the amplification process further in the nonlinear regime, we obtain that 
fluctuations of the radiation intensity and radiation power increase, and relevant probability 
distributions tend to those given by eqs. ([141) and ([T5]) . This behavior hints that the properties 
of the radiation from a SASE FEL operating in the deep nonlinear regime tend to be those 
of completely chaotic polarized light [25, 34]. For the deep nonlinear regime we find that the 
degree of transverse coherence defined by ([TOl) again tends to be an agreement with ([T6l) . 

Another practical problem refers to the probability distributions of the radiation intensity 
in the frequency domain, like that filtered by monochromator. For SASE FEL radiation 
produced in the linear regime the probability distribution radiation intensity is defined by 
gaussian statistics, and it is the negative exponential for a narrow band monochromator. 
When amplification process enters saturation regime, this property still holds for the case 




(14) 




(15) 
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of long electron pulse [18,25], and is violated significantly for the case of short electron 
bunch, of about or less than coherence length. In the latter case fluctuations of the radiation 
intensity after narrow band monochromator are significantly suppressed as it has been 
predicted theoretically and measured experimentally at the free electron laser FLASH at 
DESY operating in a femtosecond mode [26, 37]. 

All considerations presented above are related to the fundamental harmonic of the SASE 
FEL radiation. Radiation from a SASE EEL with a planar undulator has rich harmonics 
contents. Intensities of even harmonics are suppressed [49], but odd harmonics provide 
significant contribution to the total radiation power [50-54]. Comprehensive studies of 
the statistical properties of the odd harmonics have been performed in paper [45]. It has 
been found that the statistics of the high-harmonic radiation from the SASE EEL changes 
significantly with respect to the fundamental harmonic (with respect to gaussian statistics). 
Eor the fundamental harmonic the probability density function of the intensity is the negative 
exponential distribution: p(W) = {W)~^ exTp(—W/ {W)). Mechanism of the higher 
harmonic generation is equivalent to the transformation of the intensity W as z = (W)^, 
where n is the harmonic number. It has been shown in [45] that the probability distribution 
for the intensity of the n-th harmonic is given by: 

p{z) = expi-z'/y{W)) . (17) 

n{W} 

The expression for the mean value is {z) = n\{W)". Thus, the nth-harmonic radiation 
for the SASE EEL has an intensity level roughly n! times larger than the corresponding 
steady- state case, but with more shot-to- shot fluctuations compared to the fundamental [54]. 
Nontrivial behavior of the intensity of the high harmonic reflects the complicated nonlinear 
transformation of the fundamental harmonic statistics. In this case a gaussian statistics are no 
longer valid. Practically this behavior occurs only in the very end of high gain exponential 
regime when coherent radiation intensity exceeds an incoherent one. When amplification 
enters nonlinear stage, probability distributions change dramatically on a scale of the gain 
length, and already in the saturation regime (and further downstream the undulator) the 
probability distributions of the radiation intensity of higher harmonics are pretty much close 
to the negative exponential distribution [45]. 



5. Characteristics of the radiation from SASE FEL operating in the saturation regime 

In Eig.[2]we present evolution of the main characteristics of a SASE EEL along the undulator. 
If one traces evolution of the brilliance (degeneracy parameter) of the radiation along the 
undulator length, there is always the point (defined as the saturation point [34]) where the 
brilliance reaches maximum value. The best properties of the radiation in terms of transverse 
and longitudinal coherence are reached just before the saturation point, and then degrade 
significantly despite the radiation power continuing to grow with the undulator length. 

Application of similarity techniques allows us to derive universal parametric 
dependencies of the output characteristics of the radiation at the saturation point. As we 
mentioned in Section 2, within accepted approximations (optimized SASE EEL and negligibly 
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Figure 7. Degree of transverse coherence ^gat in the saturation point versus parameter e. The number of 
electrons in the coherence volume is A^c = 4 x 10^. Calculations have been performed with the simulation 
code FAST [48]. 

small energy spread in the electron beam), normalized output characteristics of a SASE 
FEL at the saturation point are functions of only two parameters: e = 27re/A and the 
number of electrons in the volume of coherence = INg\/c, where A^g = L^/\^ is 
the number of undulator periods per gain length. Characteristics of practical interest are: 
saturation length Lgat, saturation efficiency 77sat = Psat/Ph (ratio of the radiation power to the 
electron beam power Pb = 'jmc^I /e), coherence time Tc, degree of transverse coherence C, 
degeneracy parameter 5, and brilliance B^. Applications of similarity techniques to the results 
of numerical simulations of a SASE FEL [34] gives us the following result: 

4at = rL,at^2.5xe^/6xlnA^, , 

fi = P/(pPb) ^0.17/e, 

fc = puTc ~ 1.16 X a/Ih A'c X e^/^ , 

(Tu = v^/tc. (18) 

These expressions provide reasonable practical accuracy for e > 0.5. With logarithmic 
accuracy in terms of Nc characteristics of the SASE FEL expressed in a normalized form 
are functions of the only parameter e. The saturation length, FEL efficiency, and coherence 
time exhibit monotonous behavior in the parameter space of modem XFELs (e ~ 0.5 ... 5). 
Situation is a bit complicated with the degree of transverse coherence as one can see in 
Fig. Ul The degree of transverse coherence reaches a maximum value in the range of e ~ 1, 
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and drops at small and large values of e. At small values of the emittance, the degree 
of transverse coherence is limited by the interdependence of poor longitudinal coherence 
and transverse coherence [32]. Due to the start-up from shot noise, every radiation mode 
entering eq. (21) is excited within finite spectral bandwidth. This means that the radiation 
from a SASE PEL is formed by many fundamental TEMoo modes with different frequencies. 
The transverse distribution of the radiation field of the mode is also different for different 
frequencies. Smaller values of e (smaller value of the diffraction parameter) correspond to 
larger frequency bandwidths. This effect explains the decrease of the transverse coherence 
at small values of e. The degree of transverse coherence asymptotically approaches unity as 
(1 — C) oc 1/ z (X 1/ In A'c at small values of the emittance. 

In the case of large emittance the degree of transverse coherence is defined by the 
contents of higher transverse modes [34, 35]. When e increases, the diffraction parameter 
increases as well, leading to the degeneration of the radiation modes [18]. The amplification 
process in the SASE EEL passes limited number of the field gain lengths, and starting from 
some value of e, the linear stage of amplification becomes too short to provide a mode 
selection process dH). When the amplification process enters the nonlinear stage, the mode 
content of the radiation becomes richer due to independent growth of the radiation modes in 
the nonlinear medium. Thus, at large values of e the degree of transverse coherence is limited 
by poor mode selection. The degree of transverse coherence scales as ^sat (In N^/ e)^ in the 
asymptote of large emittance. To avoid complications, we present here just a fit for the degree 
of transverse coherence for the number of electrons in the coherence volume Ai'c = 4 x 10^: 

Csat ^ YToIsi^ ■ ^^^^ 

Recalculation from reduced to dimensional parameters is straightforward. Eor instance, 
saturation length is L^at — 0.6 x Lg x In N^,. Using (fTSl) and (fT9l ) we can calculate normalized 
degeneracy parameter 5 = fjC,fc and then the brilliance (fT2)) : 

[ ./""r^, , J - 4.5 X 10- X 'm^ElGeV] ^ - ^^^^ 

[sec mrad"' mm^ 0.1% bandw.) J A[A ] 

Properties of the odd harmonics of the radiation from a SASE PEL with a planar 
undulator operating in the saturation regime also possess simple features. In the case of cold 
electron beam contributions of the higher odd harmonics to the PEL power are functions of 
the only undulator parameter K [45]: 

(W^s), Kl (1^5), Kl 

Here = A'(-l)<''-"/^|J„.-i,/2(0) - Jtk+i,/2{Q)]. Q = A'V[2(1 + and h is an 
odd integer. Influence of the energy spread and emittance leads to significant decrease of the 
power of higher harmonics, up to a factor of three for the third harmonic, and a factor of up 
to ten for the fifth harmonic. Power of the higher harmonics is subjected to larger fluctuations 
than the power of the fundamental harmonic as we mentioned in the previous section. The 
coherence time at saturation scales inversely proportional to the harmonic number, while 
relative spectrum bandwidth remains constant with the harmonic number. 
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6. Estimations in the framework of the one-dimensional model 



An estimation of SASE PEL characteristics is frequently performed in the framework of the 
one-dimensional model in terms of the EEL parameter p [55]: 

1/3 

(22) 



At. 



where jo = //(27rcr^) is the beam current density, a = \f]ie^J^ is rms transverse size of 
the electron beam, and (5 is external focusing beta function. EEL parameter p relates to 
the efficiency parameter of the 3D EEL theory as = p/ B^^^. Basic characteristics of the 
SASE EEL are estimated in terms of the parameter p and number of cooperating electrons 
Nc = I /{epuj). Here we present a set of simple formulae extracted from [18,22,25]: 



The field gain length 

Saturation length 

Effective power of shot noise 
Saturation efficiency 
The power gain at saturation 

Coherence time at saturation 
Spectrum bandwidth 



-Psh 



Airp 

A,,, 



47rp 



3 + 



IniVc 
V3 



p , 

G ~ ^7V, vVhTiVe , 
1 hlnNr 



puj V 18 



(23) 



In many cases this set of formulas can help quickly estimate main parameters of a SASE EEL 
but it does not provide complete self-consistent basis for optimization of this device. 
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